2*((x - x2)*((x - x1)^2 + (y - y1)^2)*(-r12^2 + (r11 -
    ((x - x1)^2 + (y - y1)^2)^(1/2))^2 + (z - z1)^2)*(-r22^2 + (y - 
    y2)^2 + (r21 - ((x - x2)^2 + (z - z2)^2)^(1/2))^2) + (x - 
    x2)*((x - x1)^2 + (y - y1)^2)*(-r12^2 + (r11 - 
    ((x - x1)^2 + (y - y1)^2)^(1/2))^2 + (z - z1)^2)*(-r21 + 
    ((x - x2)^2 + (z - z2)^2)^(1/2))*((x - x2)^2 + (z - 
    z2)^2)^(1/2) + (x - x1)*(-r11 + 
    ((x - x1)^2 + (y - y1)^2)^(1/2))*((x - x1)^2 + (y - 
    y1)^2)^(1/2)*(-r22^2 + (y - y2)^2 + (r21 - 
    ((x - x2)^2 + (z - z2)^2)^(1/2))^2)*((x - x2)^2 + (z - 
    z2)^2) + (x - x1)*(-r12^2 + (r11 - ((x - x1)^2 + (y - y1)^2)^(1/2))^2 +
    (z - z1)^2)*(-r22^2 + (y - y2)^2 + (r21 - 
    ((x - x2)^2 + (z - z2)^2)^(1/2))^(2))*((x - x2)^2 + (z - z2)^2));


2*(((x - x1)^2 + (y - y1)^2)*(y - 
    y2)*(-r12^2 + (r11 - ((x - x1)^2 + (y - y1)^2)^(1/2))^2 +
    (z - z1)^2) + (-r11 + ((x - x1)^2 + (y - y1)^2)^(1/2))*
    ((x - x1)^2 + (y - y1)^2)^(1/2)*(y - y1)*(-r22^2 +
    (y - y2)^2 + (r21 - ((x - x2)^2 + (z - z2)^2)^(1/2))^2) +
    (y - y1)*(-r12^2 + (r11 - ((x - x1)^2 + (y - y1)^2)^(1/2))^2 +
    (z - z1)^2)*(-r22^2 + (y - y2)^2 + (r21 - 
    ((x - x2)^2 + (z - z2)^2)^(1/2))^2))*((x - x2)^2 + (z - z2)^2);

2*((x - x1)^2 + (y - y1)^2)*((z - z1)*(-r22^2 + (y - y2)^2 +
    (r21 - ((x - x2)^2 + (z - z2)^2)^(1/2))^2)*((x - x2)^2 +
    (z - z2)^2) + (-r12^2 + (r11 - ((x - x1)^2 + (y - y1)^2)^(1/2))^2 +
    (z - z1)^2)*(-r22^2 + (y - y2)^2 + (r21 - ((x - x2)^2 + (z - z2)^2)^(1/2))^2)*
    (z - z2) + (-r12^2 + (r11 - ((x - x1)^2 + (y - y1)^2)^(1/2))^2+ 
    (z - z1)^2)*(-r21 + ((x - x2)^2 + (z - z2)^2)^(1/2))*
    ((x - x2)^2 + (z - z2)^2)^(1/2)*(z - z2));